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FOREWORD 


This document presents the results of a study to develop a procedure for 
the determination of the effect of fuselage nose bluntness on the wave drag 
of supersonic cruise aircraft. The Vought Corporation Hampton Technical Center 
provides technical support to the Advanced Supersonic Technology Office, 
Aeronautical Systems Division, NASA Langley Research Center under Contract 
Number NASI -13500. The study was monitored by Mr. Vincent R. Mascitti of the 
Advanced Supersonic Technology Office. This report was prepared by 
Mr. Kenneth B. Walkley under the direction of Mr. C. W. Pearce, the Hampton 
Technical Center Advanced Aircraft Projects Supervisor. 
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A PROCEDURE FOR THE DETERMINATION OF THE 
EFFECT OF FUSELAGE NOSE BLUNTNESS ON THE 
WAVE DRAG OF SUPERSONIC CRUISE AIRCRAFT 

* By Kenneth B. Walkley 

Vought Corporation Hampton Technical Center 

SUMMARY 

The incremental wave drag penalty due to nose blunting of a fuselage has 
been investigated using a three-dimensional finite difference scheme. An 
aircraft typical of current supersonic cruise concepts has been considered. 
Computational problems with the finite difference scheme as the fuselage 
afterbody closes have been addressed. A linear theory method has been 
employed to compute the afterbody aerodynamics and effectively extends the 
finite difference scheme to closing afterbodies. Acceptable drag increments 
for various levels of nose bluntness have been demonstrated using this 
approach. 


INTRODUCTION 

The ability to predict the effect of nose bluhtness on the wave drag of 
supersonic cruise aircraft is important for the aircraft designer. Recent 
studies of sonic boom minimization (references 1 and 2) stress the importance 
of nose bluntness in reducing overpressure levels. Continuing NASA low 
speed aerodynamics tests (reference 3) are aimed at achieving low attitudes 
in the approach and landing condition for supersonic cruise aircraft. These 
low attitudes coupled with nose blunting could simplify a complex mechanical 
visor nose which' is now required for aircraft such as the Anglo-French 
Concorde. Reduced attitudes which maintain good pilot visibility could' 
result in reduced weight and complexity, shorter landing gear struts,- and 
improved performance. These potential benefits for supersonic cruise' aircraft 
will be paced, however, by the increases in wave drag associated with blunt- 
ing of the nose. 
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Many methods for the, calculation of blunt nose aerodynamics for bodies of 
revolution exist in the literature (see references 4, 5, 6, 7, and 8 for 
example-). Although the accuracy of these various methods has been substan- , 
tiated, application to specific aircraft designs at arbitrary supersonic Mach 
numbers in a straightforward manner is not usually possible. The finite 
difference relaxation method of reference 4, for example, applies to transonic 
Mach numbers only while for flight conditions below Mach 3.0, the method of 
characteristics approach of reference 5 is not generally applicable. Similarly, 
only the blunt nose of a fuselage may be analyzed using the procedures described 
in references 6 and 7. The method of reference 8 allows for an arbitrary 
afterbody shape, but no results for fuselage type bodies with closure are 
presented. All of these methods are available in computer program form, but 
both specific input requirements as well as the associated program output 
data vary considerably. Often program output is not in'a form which is 
conveniently usable by the aircraft designer. 

The three-dimensional finite difference codes of reference 9 overcome some 
of the shortcomings noted above and provide the means for a relatively straight- 
forward analysis of complex aircraft geometries including the effects of nose 
bluntness. The purpose of this paper is to present the results of an analysis 
of nose bluntness effects for a supersonic cruise vehicle using the finite 
difference codes of reference 9. An assessment of the accuracy of the codes 
is also considered through comparisons with both experimental data and a linear 
theory method. 

The author wishes to acknowledge the contributions to this document of 
Vincent R. Mascitti of the NASA Langley Research Center. Mr. Mascitti developed 
the numerical approach for -applying the Lighthill integral to arbitrary 
axisymmetric bodies. The Appendix of this paper summarizes the development 
of this method. His suggestions and guidance during the course of this effort 
are greatly appreciated. 
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TECHNICAL APPROACH 
Configuration Selection 

‘A configuration typical of current AST (Advanced Supersonic Technology) 
design concepts has been chosen as the basis for this analysis. Both the 
original sharp nosed fuselage and four blunt nosed variations have been studied 
at Mach numbers of 2.7, 2.2, and 1.7 to assess the incremental wave drag, 
penalties associated with nose .bluntness. All drag coefficients presented' 
herein are based on the configuration wing reference area and may be applied 
directly to the complete configuration . j -k ■ 

The AST-100 is a 292-passenger concept designed for Mach 2.7 cruise at 
18,288 m (60,000 ft). A complete description and-analysis of , this, configura- 
tion has been presented in reference 10, and a three-view drawing of the 
aircraft is presented in Figure 1. Of particular interest to this study is 
the fuselage. This high fineness ratio body (&/d max - 25.5) is circular in 
cross-section except in the mid-region where it is slightly elliptical. The 
fuselage is also cambered for integration with the theoretical wing root camber 
and area-ruled for minimum overall configuration wave drag at the Mach 2.7 
cruise condition. The baseline fuselage employs a sharp nose aligned with 
the freestream flow. For analysis purposes, this fuselage has been modeled 
as an uncambered body of revolution with the appropriate area distribution. 

An angle of attack of zero degrees has been assumed throughout the analysis. 

Justification of Body-Alone Analysis 

f 

The purpose of the analysis has been to evaluate the incremental drag 
penalties associated with various levels of nose bluntness relative to the 
baseline sharp nose. These penalties have been assessed through analysis of 

*■ i 

the wave drag characteristics of the .fuselage alone. The validity of this 
approach has been demonstrated experimental ly t for a wing-body in reference 11. 
At Mach 1.61 these results indicated that nose blunting had no interference 
effects on the wing at zero angle of attack. This conclusion is further 
substantiated in Figure 2 for the AST-100 where both configuration wave drag 

1 ♦ * J T , | » 

, and interference drag data are presented for a range of supersonic Mach 
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numbers. These results have been obtained using the near-field wave drag 
analysis method 'presented in reference 12. This method assumes zero-lift 
conditions and calculates the thickness and interference pressure distributions 
for the configuration with zero wing and fuselage camber. The resultant 
drag -force is obtained by integrating the pressure distributions over the 
components cross-sectional areas. Figure 2 shows that for the complete 
configuration the interference drag coefficients are bounded in magnitude by 
two drag counts ( + .'0002) whereas -for the wing-body alone, the total interfer- 
ence drag coefficient is less than one drag count (.0001) for the Mach number 
range of interest. These results indicate that much of the interference drag 
occurs between the nacelles and other aircraft components. It is assumed that 
these nacelle terms are not appreciably influenced by changes in the local 
flow at the nose, and thus a reasonable approach to the current problem is 
to consider the fuselage alone. Variations in the fuselage drag levels due 
to nose bluntness thus can be applied directly to the baseline configuration 
drag levels. 


ANALYSIS AND RESULTS 

< 

Consideration of the present need to analyze flows over blunt nosed 
bodies of revolution for which the linear theory is inapplicable has resulted 
in the choice of a three-dimensional finite difference scheme for calculation 
of the fuselage aerodynamics. The method chosen has been documented in 
reference 9 while an abbreviated description of the method and correlations 
of the theoretical results with experimental data have been presented in 
reference 13. A recent review of efforts with these finite difference codes 
as applied to supersonic cruise aircraft has been published as reference 14. 

Code Description and Capability 

The finite difference code is comprised of five separate computer programs, 
three of which have been employed herein. The QUICK code provides the means 
for modeling relatively complex geometries and providing the continuous 
analytic definition required by the other codes. Computation of the super/ 
hypersonic flows about rather arbitrary configurations over a wide range of 
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Mach numbers and angles of attack is handled by the STEIN code. Use of this 
program first requires a complete definition of the flow field at a starting 
plane location near the nose of the configuration. A routine for definition 
of this starting plane for sharp circular cones at small angles of attack 
is provided within the STEIN code itself. Starting plane data for flows over 
blunt noses is generated using the BLUNT code which is compatible with both 
the QUICK and STEIN programs (the BLUNT code is based on reference 7). Once 
the starting plane data is defined, the STEIN code computes the complete flow 
field between the configuration and the associated shock envelope from the 
starting plane to a specified end station. Additional details on the theory 
and computational techniques employed in these methods may be obtained from 
the references cited above. 

Demonstration of Code Applicability 

Blunt cone pressure distributions have been computed at Mach numbers of 
2.96, 2.30, and 1.90 using the BLUNT and STEIN codes to demonstrate their 
applicability to the present effort. The geometry of the blunted cone 
considered and the results obtained are presented in Figure 3. The experi- 
mental data have been obtained from reference 8. The agreement between the 
theoretical estimates and the measured data is quite good at all Mach numbers 
considered. Difficulty with execution of the BLUNT code occurs as the Mach 
number is decreased below 1.90, and no solution was obtained at Mach 1.50. 

Current efforts toward extending the applicability of the BLUNT code to these 
lower supersonic Mach numbers are continuing at the NASA Langley Research Center, 
but no attempt to obtain solutions for the present problem using the revised 
methods have been made. 

The capability of the finite difference code to predict afterbody aero- 
dynamics is demonstrated in Figure 4. The experimental data shown have been 
obtained by William K. Abeyounis of the NASA Langley Research Center and are 
unpublished. The agreement between the theory and experiment is excellent 
to within approximately ten percent of the body closure point. The computed 
pressure coefficients then increase rapdily, and the program actually 
terminates due to numerical errors prior to computation at the closure point. 

6 oeigw^page® 
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Afterbody closure effects in supersonic flows have been investigated by 
Meyer (reference 15) and Reyn (reference 16). Both authors point to the 
possibility of the formation of a region of subsonic flow ahead of the body 
closure point. The region is preceded by an attached tail shock and is 
followed by a sonic line. If such a subsonic region does occur on the 
fuselage under consideration, the STEIN code cannot continue to compute. 
Existence of this subsonic region thus provides a plausible explanation 
for the failure of the STEIN code near the rear, closure point. A procedure 
for effectively extending the finite difference scheme to closing afterbodies 
is presented below in the application of the theory to the AST-100 fuselage. 

Application to Present Study 

The uncambered, circular cross-section representation of the* AST-100 
fuselage has been modeled in two parts using the QUICK code as shown in 
Figure 5. A two-part definition is required to attain the desired level of 
accuracy in the model while not exceeding the QUICK code limitations on the 
maximum number of arcs allowed for the geometry definition. As Figure 5 
shows, good accuracy in the model radius distribution has been achieved using 
a combination of linear segments connected by second-order fairings. This 

i 

sharp nosed fuselage has been blunted as shown in Figure 6. The associated 
reduction in fuselage length and fineness ratio are also summarized in 
Figure 6. The body definition aft of the various blunt noses remains 
unchanged relative to the original sharp nose in all cases. 

Required starting plane solutions have been computed using the STEIN and 
BLUNT codes as previously described. The choice of the starting plane 
location is arbitrary although some care must be exercised with the blunt 
noses to insure that the axial component of the Mach number at the starting 
plane is supersonic. Pressure distributions typical of those obtained are 
shown in Figure 7. Experience with the BLUNT code has indicated that the 
axial Mach number at the starting plane must be "sufficiently" supersonic 
to achieve successful execution of the STEIN code on subsequent runs. In 
particular, it may be necessary to move the starting plane location further 
downstream for the blunt noses as the freestream Mach number is decreased in 
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order to obtain a reasonable starting plane Mach number. 

Having defined the necessary starting plane data for the various geome- 
tries, the STEIN code has been used to compute the body pressure distributions 
and associated wave drag coefficients. It should be noted that two 
executions of the STEIN code are required for computation of the flow field 
over the entire body length because of the two-segment geometry definition. 

This procedure is relatively straightforward in that the flow field at any 
intermediate body station may be saved on punched cards, magnetic tape, or 
disk, and then used to restart the code for subsequent calculations. This 
approach does increase analysis time, however, due to the relatively long 
turnaround times currently associated with the STEIN code. 

Finite Difference Code Results 

Sharp Nosed Fuselage - Computed surface pressure distributions at Mach 
numbers of 2.7, 2.2, and 1.7 are presented in Figure 8. Results obtained 
from both the STEIN code and the linear theory Lighthill integral method are 
shown. A detailed description of the Lighthill method as well as demonstrations 
of its applicability and accuracy are included in the Appendix. The STEIN 
code results have been obtained using 900 mesh points over the forward portion 
of the fuselage and 1800 over the aft section. This increase in the number 
of mesh points is required because of the relatively large distances between 
the body and shock toward the fuselage aft section. The Lighthill solution 
has been obtained using 501 points to represent the body. The QUICK geometry 
modeling has been employed for both methods. As Figure 8 shows, the pressure 
distributions computed with these two methods agree quite well over the 
majority of the body length. Both methods encounter some difficulty as the 
rear stagnation point is appraoched, however. The STEIN code result appears 
to be more sensitive to the body closure, and the computed pressure coefficient 
instability becomes more acute as the Mach number decreases. As previously 
noted, the STEIN code is unable to compute the flow field to the end of the 
body, and the program stops with a numerical error exit prior to the closure 
point. The Lighthill method does compute the pressure distribution to a 
point arbitrarily close to the body closure point, but the magnitudes of the 
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computed values can become relatively large. 

As shown in Figure 9, both the magnitude and trend of the wave drag 
coefficients computed by the STEIN code disagree markedly with results 
obtained using the Lighthi 1 1 method. This result is to be expected, however, 
as the code failure at the rear of the fuselage results in large positive 
pressures which provide an erroneous thrust. If these aft pressures (for 
fuselage stations beyond 73 m (240 ft) for example) are replaced with the 
values computed using the Lighthi 1 1 theory, results which closely agree with 
the Lighthi 11 values can be obtained as shown in Figure 9. - 

Blunt Nosed Fuselage - Figure 10 presents typical comparisons of the 
pressure distributions computed using the BLUNT and STEIN codes for the 
blunt nosed fuselage R-_--/R m - v = .246 and the original sharp nosed case. 

Note that the origin of the sharp nosed body has been shifted to the left to 
maintain geometric correspondence between the two bodies. Comparison of 
these pressure distributions shows that the blunt nosed fuselage pressure 
distributions correspond to the sharp nosed solution at the same Mach number 
for points downstream of the nose region, and that the effects of nose blunt- 
ing tend to be localized on the forebody. This result further confirms the 
previous conclusion that changes to the nose geometry do not result in 
appreciable interference on the wing. 

The total drag for the blunt nosed fuselages has been computed using the 
BLUNT and STEIN codes pressures aft to approximately station 73* m (240 ft) 
and then substituting the Lighthi 11 values as was done for the sharp nosed 
cases. The results of this analysis are summarized in Figure 11 where the 
incremental wave drag coefficient due to blunting is presented. There 
appears to be a negligible drag penalty for noses with bluntness values less 
than about 0.15, but the incremental drag then increases at a fairly high 
rate. 


The relative magnitudes of these drag increments are illustrated in 
Figure 12 where typical values of skin friction, wave drag, and drag-due-to- 
lift are presented for the AST-100 configuration. Reasonable nose bluntness 
levels would results in additional drag penalties of less than ten percent 
of the configuration total drag at the Mach 2.7 cruise condition. 
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CONCLUDING REMARKS 


The incrementa-1 wave drag penalty due to nose blunting of a fuselage has 
been investigated using a three-dimensional finite difference scheme. An 
aircraft typical of current supersonic cruise concepts has been considered. 

The analysis assumes that the incremental drag value may be determined by 
considering the fuselage alone. 

The finite difference scheme provides an exact inviscid solution for the 
blunt nose and forebody, but computational problems occur near the afterbody 
closure point. The existence of a region of subsonic flow has been postulated 
as a possible cause of the failure of the finite difference scheme in this 
area. A linear theory method has been used to compute the afterbody aero- 
dynamics for the sharp nosed fuselage. These results may then be matched 
with the finite difference scheme ahead of the afterbody to obtain the complete 
fuselage pressure distribution and wave drag. This approach is also valid for 
the blunt nosed fuselages since the effects of blunting are shown to be 
localized on the forebody. 

Acceptable drag increments for various levels of nose bluntness have been 
demonstrated for a typical supersonic cruise concept using these modified 
finite difference code pressure distributions. The relative importance of 
these drag increments will generally depend on the specific vehicle applica- 
tion, however. 

A possible extension of the present method would utilize the finite 
difference codes to compute the aerodynamics for the blunt nose with the linear 
theory method providing pressures for the remainder of the body. This approach 
would probably be the most convenient and computationally efficient. 
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NOTE DIMENSIONS SHOWN IN METERS WITH FEET 
IN PARENTHESIS EXCEPT AS NOTED 
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Figure 1. The AST-100 Configuration 
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Figure 2. Near- Field Wave Drag for the 1 AST-100, 
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Figure 3. Accuracy of the Finite Difference Codes as 
Applied to a Typical Blunted Cone 
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Figure 4. Afterbody Aerodynamics Using the 
Finite Difference Code 
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Figure 5. Quick Code Modeling of the AST-100 
Equivalent Fuselage 


RADIUS 


AST-100 Equivalent Fuselage 


(ft) M 



r nose 

Fuselage 
Length 
m (ft) 

‘Fineness 

Ratio 

Sharp 

96.0 

(315.0) 

'25. 6 

.137 

94.7 

(310.8) 

25.1 

.246 • 

93.7 

',(307.5) , 

24 ..9 , 

.359 

92.i ; 
'(302.3) 

24.5 . 

..'460 

90.5 

(296.8) 

24.-0 ; 


Figure .6. Blunt Nose Geometry Definitions 
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Figure 7. Typical Blunt Nose Pressure Distributions Obtained 
from the Finite Difference Code 
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Figure 8. Comparison of Computed Pressure Distributions 
fop the Sharp Nosed Fuselage 
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Figure 8. Continued 
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Figure 9. Wave Drag Variation with Mach Number 
for the Sharp Nosed Fuselage 
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APPENDIX 


Lighthill Integral Method for Axi symmetric Bodies 


A method developed by M. J. Lighthill has been used to calculate body 
pressure distributions and the associated zero-lift wave drag coefficients 
for a typical sharp-nosed AST (Advanced Supersonic Technology) fuselage. The 
purpose of this appendix is to discuss the numerical approach to the solution 
of the Lighthill equation, to validate the method through correlation of results 
obtained for several minimum drag bodies of revolution with appropriate 
experimental data, and to discuss the convergence characteristics of the 
method. A listing of both the FORTRAN computer program written to implement 
the technique and output for a typical case are also included. 


Lighthill Equation and Numerical Approach 

The equation for the surface pressure coefficient on a body of revolution 
has been shown by Lighthill in reference A-l to be 

x 


C.-1 

P IT 


f UJ &- - 


(A-l) 


where 

x = body field station 
U(Z) = decay function 

X — t 

Z = position function, = 

t = x variable of integration 

6 = Mach number parameter, = - 1 

R(t) = body radius at t 

S'(t) = first derivative of body cross-sectional area S at t 
R‘(x) = first derivative of body radius R at x 

Equation (A-l) is valid for any slender body of revolution whether it is 
smooth or not. It is easily evaluated numerically because the integrand is 
.without singularities. The decay function U(Z) is shown in Figure (A-l) where 
the function 1/z is also presented for comparison. U(Z) is zero for z < 0, 
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unity at z: = 0, and- falls asymptotically to zero as 1 Z ->■ The, d'eveilopment: 
of U(Z.) has- been' presented- by LighthiTl in reference A-2: while, reference- A--? 
further discusses, the implementation in the present problem.. 

Ilf a given, body of revolution! is divided into* a finite. 1 number of stations: 
CX,R(X)],- equation (A-T) may be rewritten as:' 

Cp-f^ D k (S ’k - S'k-T )5 - tRi "l )!2 ’ ^ 


where the integral has been, replaced by a summation; from: t'he; first body^ station* 
K = 1 to the' Fth point on- the body. Similarly',. dSZ'ft)} has' Been 1 written* as; 
the difference S :1 ^ - S.'^ y. A, mean’ value of the. decay, function- onta- givem 
x- interval may be. obtained using. 1 ogarithms as follows;:; 


Iog-O k 4. 


i,ns 


k.-T; 


109 
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or. 


1 

2 

0, = e L 


09 lifer 


(M); 


thus. 


U k = 



m_) 


\6Rk l 

\ P R k-l/ 


(A-5:)) 


The use of logarithms in th.is development, is; advantageous; because, the’ 
function log [U(*Z)/eR(.t),] i, s approximately linear in. Z... 

The expression 1 for numerically evaluating the- pressure coefficient at a; 
given point I on a slender body of revolution is thus 1 :: 


v-f < s 'k - s w - ( r, p 2 


(A-6.J5 


sssss* 
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where the last term may be evaluated simply as: 


R' 


I 


= Rj - _ R I .-1 
XI - Xi_i 


(A— 7 ) 


Note in equation (A-6) that the radius R must always be non-zero. The 
solution for a given closed body must start just aft of the nose and terminate 
just ahead of the afterbody closure point. Let the solution begin at X=X-j 
such that the first two body stations Xq and X-j define an initial cone of 
half-angle ^ . The contribution of the initial cone at X=X-j may be evaluated 
as follows: At X Q the radius Rq = 0 and 


S 


0 


2nR r 



(A— 8 ) 


Now as t Xq the radius R(t) 0 and Z «>. Thus, U(Z) 1/Z: 


_ 1 _ bRq = bRq. 
0 Z X ] -t X-| 


(A— 9) 


or, 


Mil- 1_ 
3Rq X 1 


as t ■+■ X p (X-j - t)/BR(t) ■+■ 0 and 


= 1 


CA-10) 


(A-ll) 


such that 


Also, 


Ul _ 1 

pR-j pR-j 


(A— 1 2 ) 





= 2ttRi 6 -j 


(A-T3) 
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Substituting these relations into equation (A-6) gives 




( 2ir 6 -j 2 Xi ) - 6j 2 


(A— 1 4) 


or. 





(A-15) 


where Cp^ is the pressure 
result compares favorably 


coefficient at x = x-j due to the conical nose, 
with the slender body result (reference A-3) 


This 


C 


P 



2 l0 9 h - 1 


(A-16) 


No special treatment is required if the afterbody closes. The calculation is 
stopped just forward of the closure point such that the neglected drag contri- 
bution is insignificant. 


Once the body pressure distribution is known, the drag coefficient is 
calculated as 


C D 2S 


REF 


k=N 


k=l 


(R k 2 - R k-1 } 


(C p + Cp ) 
K k K k-1 


(A— 1 7 ) 


where is the reference area -and N is the total number of stations used 
to represent the body. Note again that insignificant drag contributions 
associated with the initial cone and afterbody closure point are omitted. 

This technique for calculating the pressure distribution and wave drag 
coefficient of a slender body of revolution has been programmed for the CDC 
Cyber series computers using the FORTRAN Extended language. A listing of this 
code and output for a typical case are presented at the end of this appendix. 
It should be noted that the execution times for this method are very fast (on 
the order of a few seconds) and that the solution is inherently free of 
potential numerical instabilities. 
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Accuracy and Convergence Characteristics 


Typical convergence characteristics of the Lighthill method are shown in 
Figure A-2. The smooth body of revolution considered has a fineness ratio of 
8.0 and the freestream Mach number is 3.0. As the figure indicates, the 
Lighthill integral converges rapidly to the final solution and agrees quite 
well with the method of characteristics solution. 

Three minimum drag bodies of fineness ratio 7, 10, and 13 have been studied 
using the Lighthill method. As shown in Figures A-3 and A-4, the predicted 
pressure distributions and wave drag coefficients show excellent agreement with 
both the method of characteristics solution and the experimental data (reference 
A-4) for a wide range of Mach numbers. 

Results obtained for a typical closing afterbody are presented in Figure 
A-5. The experimental data have been obtained by William K. Abeyounis of the 
NASA Langley Research Center and are unpublished. Although the Lighthill 
theory overpredicts the afterbody pressures, the agreement is still reasonable. 
Note in particular that a finite pressure coefficient is obtained at the 
point where the solution is terminated (x/£ = 0.997). 

The results presented above have indicated both good accuracy and conver- 
gence characteristics of the method when applied to smooth analytic shapes. 
Analysis of typical area-ruled fuselages has pointed to some lack of convergence 
in the computed wave drag coefficient, however. As shown in Figure A-6, 
satisfactory convergence has not been achieved using as many as 501 points 
to define the body radius dsitribution. The extrapolation technique 
illustrated in the figure can be used to estimate the converged results as 
shown. The resulting wave drag variation with Mach number is also shown in 
Figure A-6. 

In conclusion, the Lighthill technique presents a rapid method for 
determination of pressure distributions and wave drag coefficients for bodies 
of revolution in supersonic flow. Care must be taken, however, to insure 
that converged results are obtained. 
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Figure A-5. Correlation of Afterbody Pressure 
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APPENDIX 


PROGRAM DR AG (INPUT»0UTPUT#TAPE5*INPUT»TAPE6» OUTPUT) 

CALCULATES PRESSURE DISTRIBUTION AND WAVE DRAG FOR A SLENDER 
POINTED BODY OF REVOLUTION USING THE LIGHTHILl INTEGRAL FOR 
NON-SMOOTH BODIES 


********************* INPUT REQUIREMENTS *************************** 

CARD 1 TITLE INFORMATION - FORMAT! 8A10) 

CARD 2 NMACH* NPTS - F0RMAT(2I3) 

NMACH-NUMBER OF MACH NUMBERS ULE.10) 

NPTS -NUMBER OF POINTS DEFINING BODY GEOMETRY 

CARD 3 AM - ARRAY OF MACH NUMBERS - F OR M AT ( 1 0F7 .0 ) 


NOTE* SUBROUTINE GEOM PROVIDES ( STA TI ON , RAD IUS ) 

DATA AND MAY BE SET UP TO COMPUTE OR 
READ AS REQUIRED FOR A GTVEN PROBLEM. 

*********************************************************************** 
DIMENSION X(5Q1 )* R ! 501 )* S t 501 )*RP (501), SP I 501)* CP <501) 

DIMENSION ABC ( 8 ) # AM ( 10 ) 

READ TITLE 

500 READ( 5*10) ABC 
10 FORMAT ( 8 A10 ) 

IF ( EOF ( 5 ) ) 20# 30 
20 STOP 
30 CONTINUE 

READ NMACH# NPTS 

READ( 5*40) NMACH. NPTS 
40 FORMAT (213) 

READ ARRAY OF MACH NUMBERS 


READ! 5*50) < AM { I ) , 1-1# NMACH ) 

50 FORMAT !1GF7.0) 

GENERATE OR READ BODY GEOMETRY 

CALL GEOM!X*R#S#S REF# NPTS) 

CALCULATE BODY RADIUS AND AREA DERIVATIVES - OMIT LAST BODY 
STATION IF RINPTS)»0. 

IF(R(NPTS).EQ.O.) NPTS-NPTS-1 

ORIGINAL PAGE IB 
OP POOR QUALITY! 
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PI-3.1A1592'654 
RP{ 1) »RH)7XU) 

$P(T) ■2«*PI\*R( 1 )*RP {II 
C. 

DQ 60 K»2* NPTS- 
KK-K-l 

DELX-X(K)-XCKK) 

RP(K),«(R(K)-R(KK) )/DELX 
SP(K)»2.*PI*R(K)*RP(.K) 

60 CONTINUE. 

MACK NUMBER, LOOP FOR GIVEN' <X,R)’ GEOMETRY 

00 1U00 L-lrNMACH. 

BETA* SORT { AM ft ) *AM ( L ) -1. ) 

CPVAC*-2./(I.4*AM(L>*AMlLn 

CALCULATE PRESSURE DISTRIBUTION AND WAVE DRAG 
CALL CPDRAG(8ETA» X# Rj-R'P*>S,P> NPTS» SREF» C'P* CD-OV 
WRITE OUTPUT FOR CURRENT M‘ACH NUMBER 

CALL OU<TLX>R*S»RP,SP'j CP,NPTS,COO>SREe* A'M(L), ABC*CPVAC) 
C 

lUOO CONTINUE 
C' 

C NEXT GEOMETRY 

C 

GO TO 500 
C 

END 
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SUBROUTINE GEOM { X# R# S> SREF# NPTS ) 

TO PROVIDE THE REQUIRED GEOMETRY DATA - THIS ROUTINE 
MUST EITHER CALCULAT£*READ#OR EXPLICITLY DEFINE THE 
t ST AT I ON# RADIUS) DATA AT NPTS POINTS 

X(l) AND R C 1 ) MUST BE GREATER THAN ZERO 

DIMENSION X(NPTS )#R(NPTS)#rS{ NPTS J 

GENERATE HAACK-ADAMS BODY IL/DMAX«13) NASA TN D-3163 

PI-3.1A159265A 

RMAX-1.385 

XL*36« 

XINC* • 005 
XU)«.001 
X(2)**005 

DO 10 I» 3# NPTS 
X(I )»X ( I-U + XINC 
10 CONTINUE 

DO 20 I»1#NPT5 
XDEL«2.*XCI)-1* 

T1-.707*U.-XDEL*XDEL)**1.5 

T2».1693A*XDEL*SQRT(1.-XDEL*XDEL) 

T3«.1693A*AC0S(-XDEL) 

R(I)«SQRT(T1+T2+T3)*RMAX 

SU)-PI*R(I)*R(I) 

X(I)«X(I)*XL 

20 CONTINUE 

DEFINE REFERENCE AREA 

SREF-PI*RMAX*RMAX 

RETURN 

END 
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SUBROUTINE C PDRAG ( BETA, X, R, RP, SP, NPTS, SREE,C P, C DO) 

TO COMPUTE THE BODY PRESSURE DISTRIBUTION AND 
WAVE DRAG COEFFICIENT USING THE LIGHTHILL INTEGRAL 

DIMENSION X.(NPT.S).,R(NPT-S),-RP(NPTS),SP(NPTS),CP(-NPTS-) 

DIMENSION PSIT(36),UTAB(36) 

DATA PSIT/O., .2, .4, .6, .8,1. ,1.2, 1.4, 1.6, 1.8, 2., 2. 2, 2. 4, 2 .6, 2*6, 3.0 
*,3.2, 3. 4, 3. 6, 3. 8, 4. ,4. 4, 4. 3, 5.2, 5. 6, 6.,6.4»6.9,*7.2,7.6>8«,8»4,8.8, 
♦9*2, 9.6,10./ 

DATA UTAB/1.,. 90703,. 82646, .75621,. 694 62,. 6403,4,. 59229,. 54960# 

*. 51149,. 47737, .44672, .*1907, . 39408, . 37140, .3508 0, .33201, .3148 3, 

*.29909, .28 464, .27134, .25906, .23721,. 2 1840, .20209, .18 785, .17534, 

*.16428, .1544 5, .14567, .13778, • 13068, • 1 2424, . 11 839, .11304, • 10815, 
4,10366/ 

PI-3.141592654 

PRESSURE COEFFICIENT AT FIRST BODY STATION 
DUE TO INITIAL CONE 

CP(1)-RP(1)*RP( 1) 4(2. /SORT ( BETA*R P(1 ) ) -1. ) 

DO 10 1-2, NPTS 
ULAST- 1 • /X ( I ) 

SUM2-0. 

DO 20 K- 1, I 

PSI-(X(I)-XCK))/(BETA4R(K) ) 

IF< PSI.GT.10. ) GO TO 30 

CALL FTLUP(P$I,U,1,36,PSIT,UTAB) 

GO TO 40 
30 U-l./PSI 
40 CONTINUE 

IFCK.EQ.l) SS-O. 

IF(K.GT.l) SS-SP(K-l) 

SUM2»SUM2+SQRT{ULAST4U/(BETA4R(K)))*{SPCK)-SS)/PI 

ULAST-U/(8ETA4R(K>) 

20 CONTINUE 

CP(I)-SUM2"'RPCI)*RPCI) 


10 CONTINUE 

COMPUTE WAVE DRAG 
D0Q«CPU>4PI4RU)*R{1) 

C 

DO 50 1-2, NPTS 

D0Q«D0Q+Pi4(Rm4Rm-Rci-i)*RU-in/2.4(CP{n+cp<i-i) > 

50 CONTINUE < 

C 

CDQ-DOQ/SREF 

C 

RETURN 

END 
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SUBROUTINE OUT C X, R# S, R P# $ P#C P# NPTS# COO# $REF# XM , ABC#CPVAC> 

TO WRITE GEOMETRY AND AERODYNAMIC DATA 

DIMENSION X ( NPTS ) #R ( NPTS )# S (NPTS )# RP (NPTS ),SP( NPTS)# CP (NPTS) 
DIMENSION ABC ( 8 ) 

WRITE (6/10) ABC# XM# CPVAC 
10 FORM AT (1H1#1X#8A10#//#13X,#7H MACH * »F10.7# 

*10X> 12H VACUUM CP * #F10.7» 

*//#8X#lHX»15X#lHR#15X# 

*1HS#15X#5HDR/DX, II X# 5HDS /DX# 11X# 2HC? » / / ) 

DO 20 1*1# NPTS 

WRITE <6# 30) X(I)#R(I)#S(I)»RP(I),SPm#CP(I) 

30 F0RMAT(6F16.8) 

20 CONTINUE 

WRITE (6# 40 ) XM#CDO# SREF 

40 FORMATt//# & X# 11HMACH NO. * , F10. 6#5X# 10HCD WAVE - » F10. 6# 5X# 7HSREF 
* * # F10. 5 ) 

RETURN 

END 


SUBROUTINE FTLUP ( X# Y# M*N# VAR I # VARD) 

********* DOCUMENT DATE 7/7/69 SUBROUTINE REVISED 7/7/69 *♦**♦♦*** 

* MODIFICATION OF LIBRARY INTERPOLATION SUBROUTINE FTLUP 
DIMENSION VARI(N)#VARD(N)#Vm»YY<2> 

DIMENSION 11(43) 

A 

* INITIALIZE ALL INTERVAL POINTERS TO -1.0 FOR MONOTONICITY CHECK 
DATA (II(J)#J*l»43)/43*-l/ 

MA* I AB S ( M ) 

* 

* ASSIGN INTERVAL POINTER FOR GIVEN VARI TABLE 

* THE SAME POINTER WILL BE USED ON A GIVEN VARI T ABLE EVERY TIME 
L I* MOD (LOCF(VARKl) )#43) + l 

I* II ( LI ) 

IF (I.GE.O) GO TO 10 
IF (N.LT.2) GO TO 10 


* MONOTONICITY CHECK 

IF (VARI(2)-VARI(1)J 1#1#3 

* ERROR IN MONOTONICITY 

2 K*LOCF ( VARI ( 1 ) ) 


ORIGINAL PAGE IS 
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PRINT 102# J#K#( VARI ( J') # J *1# N ) # ( VARD ( J > » J*1 »N) 


102 FORMAT ( 1H1# * TABLE BELOW OUT OF ORDER FOR FTLUP AT POSITION 


1# 15# /* X TABLE IS STORED IN LOCATION *# 06# // ( 8G 15. 8 ) ) 


STOP 


* 
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* MONQTQNIC DECREASING 

i DO 5 J«2#N 

IF {VARI(J)-VARI(J-1))5,2#2 

5 CONTINUE 

GO TO 10 _ 

'* MONOTONIC" INCREASING 
3 DO 6 J-2»N 

IF ( VARI(0>“VARI( J-l) )2*2#6 

6 CONTINUE 

* 

+■ INTERPOLATION 

10 IF (I.LE.O) I»1 
IF (I.GE.N) I*N— 1 
IF (N.LE.l) GO TO 8 
IF (MA.NE.O) GO TO 99 

* ZERO' ORDER 

8 Y* VARD { 1 ) 

GO TO 800 

* LOCATE I INTERVAL ( X ( I > . LE . X. LT. X ( 1 + 1 > » 

99 IF ((VARI(I)-X>+(VARI(I+1)-X>) 61#61>40 

* IN GIVES DIRECTION FOR SEARCH OF INTERVALS 

AO IN«SIGN(1.0#(VARI{I+l)-VARim)*CX-VARim )) 

* IF X OUTSIDE ENDPOINTS# EXTRAPOLATE FROM END INTERVAL 
41 IF ( ( I + IN) .LE.O) GO TO 61 

IF m+IN).GE.N) GO TO 61 
I *I+IN 

IF ( (VARICI)-X)*(VARI(I+1)-Xn 61#61#A1 
61 IF (MA.EQ.2) GO TO 200 

* 

* FIRST ORDER 

Y»<VARDII)*IVARI{I+l)-X)«VARD(I + l)*(VARim-X,) ) / (VARIU + 1 J-VARI (I) 
1 ) 

GO TO 800 

* 

* SECOND ORDER 

200 IF (N.EQ.2) GO TO 2 

IF n.EQ.(N-l)) GO TO 209 
IF (I.EQ.l) GO TO 201 

* PICK THIRD POINT 

SK- VARII I+l )-VARI ( I ) 

IF <<SK*<X-VARI<I-1) ) ).LT.<SK*(VARItT + 2)-XJ) ) GO TO 209 

201 L«I 

GO TO 702 
209 L«I-1 

702 VU)-VARI(U-X 
V ( 2) ■ VARI < L + l) “X 
V ( 3 )* VAR 1 1 L+2 )~X 

YYC1)»(VAR!DCL)*VC2)-VARD(L+1)*V<1J 1/fVARHL + l >-VARI(L) ) 

, YT(^2)«(VARD.(L+1)*VC3)-VARD(L+2)*V(2) »/ (VARKL + 2 )-VARI(L + l>> 

Y* ( YY< 1 J + V ( J 3)-YY( 2 ) *V ! 1 ) ) / ( VAR I ( L+2) -V API ( L) ) 

800 II( LI ) *1 
RETURN 
END 
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HAACK-AOAMS BODY (NASA TN D-3163) L/DMAX-13 

MACH • 2.5300000 VACUUM CP 


X 

R 

S 

.03600000 

.01993350 

.00124829 

.18000000 

.06647351 

•01388164 

•36000000 

.11141983 

.03900092 

•54000000 

•15051021 

.07116751 

.72000000 

.18612239 

.10882962 

.90000000 

.21928216 

.15106243 

1.08000000 

.25055605 

.19722396 

1 .26000000 

.26030064 

•24683006 

1.44000000 

.30875963 

♦29949592 

1.62000000 

.33610943 

.35490431 

1.80000000 

•36246320 

.41278664 

1.98000000 

.38798475 

.47291074 

2.16000000 

.41269716 

.53507275 

2.34000000 

*43668830 

.59909126 

2.52000000 

.46001457 

•66460311 

2.70000000 

.48272354 

.73206027 

2.88000000 

.50485586 

.80072739 

3.06000000 

•52644660 

.87067992 

3.24900000 

.54752633 

.94160261 

3.42000000 

.56812190 

1.01398828 

3.60000000 

.58825709 

1.08713684 

3.78000000 

.60795308 

1.16115448 

3.96000000 

.62722886 

1.23595296 

4.14000000 

.64610153 

1.31144903 

4.32000000 

•66458658 

1.38756395 

4.50000000 

•68269809 

1.46422307 

4.68000000 

.70044891 

1.54135546 

4.86000000 

.71785084 

1.61889356 

5.04000000 

.73491469 

1.69677296 

5.22000000 

.75165048 

1.77493212 

5.40000000 

.76806743 

1.85331214 

5.56000000 

.78417412 

1.93185661 

5.76000000 

.79997854 

2.01051142 

5.9 4QOOC0Q 

.81548810 

2.08922458 

6.12000000 

.83070975 

2.16794615 

6.30000000 

.84564999 

2.24662803 


.2285714 


DR/DX 


.55370823 

.32319454 

.24970176 

.21716877 

.19784548 

.18422994 

.17374384 

.16524769 

.15810554 

.15194334 

.14652090 

.14167529 

.13729118 

.13328409 

.12959038 

.12616099 

.12295733 

.11994855 

.11710959 

.11441984 

.11186217 

.10942218 

.10708766 

.10484817 

.10269471 

.10061949 

.09861569 

.39667736 

.09479921 

.09297657 

.09120529 

.08948164 

.08780229 

.08616423 

.08456474 

.08300134 


OS/OX 


.06934966 

.13496717 

.17480907 

.20537292 

.23136871 

.25381766 

•27352320 

.29103107 

•30672360 

.32087969 

.33370853 

.34537320 

.35600326 

.36570403 

.37456244 

.38265152 

.39003332 

.39676120 

.40288150 

.40843479 

.41345693 

.41797981 

.42203196 

.42563907 

.42832442 

.43160917 

.43401263 

.43605257 

.43774531 

.43910595 

.44014851 

.44088601 

.44133062 

.44149370 

.44138592 

.44101731 


CP 


.23779907 

.17448230 

.12797003 

.10518836 

.09322190 

.08463309 

.07795337 

.07254449 

.06798630 

.06405892 

.06025137 

.05727428 

.05456474 

.05209337 

.04983358 

.04772650 

.04577377 

.04380332 

.04210186 

•04048905 

.03895726 

.03750548 

.03612562 

.03471415 

•03345510 

.03224592 

.03108512 

.02996475 

.02881657 

.02777506 

.02676532 

.02578970 

.02463875 

.02386862 

.02297403 

.02210170 
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X 

R 

S 

DR/DX 

DS/DX 

CP 

6 • 48000000 

.86031492 

2.32522391 

.08147181 

.44039733 

.02125508 

6.66000000 

.87471025 

2.40368917 

.07997408 

.43953489 

.02042653 

6.84000000 

.88884139 

2.48198075 

• 07850630 

.43843844 

.01957971 

7.02000000 

.90271340 

2.56005709 

.07706673 

.43711600 

.01879361 

7.20000000 

.91633109 

2.63787804 

.07565381 

.43557517 

.01802309 

7. 38000000 

.92969898 

2.71540481 

.07426609 

.43382318 

•01726964 

7.56000000 

.94282138 

2.79259988 

.07290221 

.43186695 

.01653169 

7. 74000000 

.95570235 

2.86942698 

.07156094 

.42971304 

.01578080 

7.92000000 

.96834575 

2.94585098 

.07024112 

.42736777 

.01507318 

8.10000000 

.98075525 

3.02183787 

.06894170 

.42483716 

.01438148 

8.28000000 

.99293436 

3.09735474 

.06766167 

.42212700 

.01370177 

8.46000000 

1.00468638 

3.17236967 

.066400U 

.41924283 

.01303404 

8.64000000 

1.01661449 

3.24685175 

.06515616 

.41619001 

.01235650 

6 , 82000000 

1.02812171 

3.32077101 

.06392901 

.41297365 

.01171281 

9.00000000 

1.03941093 

3.39409840 

.06271789 

.40959873 

.01108255 

9.16000000 

1.05048491 

3.46680574 

.06152211 

.40607001 

.01046145 

9.36000000 

1.06134629 

3.53886570 

.06034100 

•40239211 

•00983003 

9.54000000 

1.07199759 

3.61025180 

.05917393 

.39856951 

•00923143 

9.72000000 

1.08244125 

3*60093832 

.05602031 

.39460651 

•00864058 

9. 90000000 

1.09267958 

3.75090032 

.05687959 

.39050731 

•00806040 

10. 08000000 

1.10271480 

3.82011362 

.05575125 

.38627597 

.00748769 

10.26000000 

1.11254906 

3.88855474 

.05463479 

.36191642 

.00690721 

10.44000000 

1.12218442 

3.95620093 

.05352975 

.37743252 

•00635313 

10.62000000 

1.13162284 

4.02303010 

.05243568 

.37282797 

•00580667 

10.80000000 

1.14086623 

4.08902084 

.05135217 

.36810642 

.00526767 

10.98000000 

1.14991642 

4.15415237 

.05027882 

.36327139 

.00472271 

11.16000000 

1.15877516 

4.21840457 

.04921525 

.35832634 

.00419875 

11. 34000000 

1.16744416 

4.28175791 

.04816111 

.35327462 

•00368287 

11.52000000 

1.17592505 

4.34419346 

.04711605 

.34811953 

.00317380 

11*70000000 

1.18421941 

4.40569290 

.04607976 

.34286427 

.00265873 

11.88000000 

1.19232875 

4.46623846 

.04505192 

.33751199 

•00216380 

12.06000000 

1.20025456 

4.52581295 

.04403225 

.33206576 

.00167461 

12.24000000 

1.20799824 

4.58439971 

•04302046 

.32652859 

.00119136 

12.42000000 

1.21556117 

4.64198264 

.04201629 

.32090344 

.00070445 

12.60000000 

1.22294468 

4.69854615 

.04101948 

•31519320 

.00023429 

12.78000000 

1.23015004 

4.75407517 

.04002980 

.30940072 

-.00023048 

12.96000000 

1.23717850 

4.80855516 

. 03904700 

.30352879 

-.00068892 

13.14000000 

1.24403126 

4.86197206 

.03807087 

.29756014 

-.00115166 

13.32000000 

1.25070947 

4.91431230 

.03710119 

.29155749 

-.00159905 

13.50000000 

1.25721427 

4.96556280 

.03613777 

.28546348 

-.00204071 

13.68000000 

1.26354675 

5.01571095 

.03518041 

.27930074 

-.00247773 

13.86000000 

1.26970795 

5.06474462 

.03422892 

.27307163 

-.00291756 

14.04000000 

1.27569891 

5.11265214 

.03328312 

.26677931 

-.00334344 
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-21 ORIGINAL PAGE IS 

OP POOR QUALITY 


X 

R 

5 

DR/DX 

DS/DX 

CP 

14.22000000 

1.28152063 

5.15942229 

.03234285 

♦26042566 

-.00376470 

14.40000000 

1.28717406 

5.20504428 

.03140794 

.25401336 

-.00418118 

14.56000000 

1.29266014 

5.24950781 

.03047823 

.24754489 

-.00459952 

14.76000000 

1.29797978 

5.29260298 

.02955357 

.24102262 

-.00500601 

14.94000000 

1.30313387 

5.33492034 

.02863383 

.23444895 

-.00540770 

15.12000000 

1.30812326 

5.37585065 

.02771885 

.22782624 

-.00580471 

15.30000000 

1.31294880 

5.41558593 

.02680851 

,22115682 

-.00620333 

15.48000000 

1.31761128 

5.45411737 

.02590269 

.21444301 

-.00659107 

15.66000000 

1.32211150 

5.49143743 

.02500124 

.20768710 

-.00697401 

15.84000000 

1.32645023 

5.52753874 

.02410407 

.20039137 

-.00735828 

16.C200000J 

1.33062823 

5.56241435 

.02321107 

.19405806 

-.00773264 

16.20000000 

1.33464621 

5.59605773 

.02232211 

.18718942 

-.00810216 

16.38000000 

1.33850489 

5.62846273 

.021437U 

.18028765 

-.00846753 

16.56000000 

1.34220496 

5.65962361 

.02055595 

.17335498 

-.00883372 

16.74000000 

1.34574710 

5.68953504 

.01967856 

.16639359 

-.00919060 

16.92000000 

1.34913197 

5.71819206 

.01880483 

.15940565 

-.00954326 

17.10000000 

1.35236021 

5.74559013 

.01793468 

.15239334 

. -.00989612 

17.28000000 

1.35543246 

5.77172506 

.01706804 

.14535882 

-.01024053 

17.46000000 

1,35834932 

5.79659310 

.01620482 

.13830424 

-.01058068 

17.64000000 

1.36111141 

5.82019084 

.01534494 

.13123173 

-.01091665 

17.82000000 

1.36371932 

5.84251529 

.01448835 

.12414343 

-.01125225 

18. GOOQOOOO 

1.36617361 

5.86356383 

.01363496 

♦11704147 

-.01157985 

18.18000000 

1.36847466 

5.88333423 

.01278473 

.10992798 

-.01190341 

18.36000000 

1.37062363 

5.90182464 

.01193758 

.10280507 

-.01222596 

18.54000000 

1.37262045 

5.91903359 

.01109347 

.09567486 

-.01254112 

18.72000000 

1.37446587 

5.93495999 

.01025234 

.08853948 

-.01205234 

16.90000000 

1.37616042 

5.94960316 

.00941415 

.08140103 

-.01315940 

19.08000000 

1.37770461 

5.96296276 

.00857884 

.07426163 

-.01346522 

19.26000000 

1.37909896 

5.97503886 

.00774638 

.06712340 

-.01376427 

19.44000000 

1.38034397 

5.98583192 

.00691673 

.05998846 

-.01405903 

19. 62000000 

1.38144014 

5.99534275 

.00608985 

.05285894 

-.01435216 

19.80000000 

1.38238797 

6.00357258 

.00526572 

.04573696 

-.01463071 

19.98000000 

1.38318795 

6.01052301 

.00444430 

.03862467 

-.01492108 

20.16000000 

1.38384055 

6.01619603 

.00362558 

.03152419 

-.01520115 

20.34000000 

1.38434627 

6.02059401 

.00280954 

.02443769 

-.01547506 

20.52000000 

1.38470558 

6.02371972 

.00199616 

.01736732 

-.01574472 

20. 7C000000 

1.38491896 

6.02557632 

.00118543 

.01031525 

-.01601000 

20.88000000 

1.38498688 

6.02616737 

♦00037734 

.00328366 

-.01627290 

21.06000000 

1.38490982 

6.02549680 

-.00042811 

-.00372525 

-.01652971 

21.24000000 

1.38468825 

6.02356898 

-.00123092 

-.01070927 

-.01678219 

21.42000000 

1.38432266 

6.02038864 

-.00203107 

-.01766619 

-.01703163 

21.60000000 

1.38381352 

6.01596095 

-.00262857 

-.02459376 

-.01727528 

21.78000000 

1.38316130 

6.01029147 

-.00362340 

-.03148973 

-.01751424 
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X 

21 * 96000000 
22 *14000000 
22. 32000000 
22*50000000 
22.68000000 
22.86000000 
23. 04000000 
23.22000000 
23.40000000 
23.58000000 
23.76000000 
23.94000000 
24.12000000 
24.30000000 
24.48000000 
24.66000000 

24.64000000 
25.02000000 
25.20000000 
25.38000000 
25.56000000 
25.74000000 
25.92000000 
26.10000000 
26.28000000 
26.46000000 

26.64000000 
26.82000000 
27.00000000 
27.18000000 
27.36000000 
27.54000000 
27.72000000 
27.90000000 
26.08000000 
26.26000000 
26.44000000 
28.62000000 
26.80000000 
28.98000000 
29.16000000 
29.34000000 
29.52000000 


R 

1.38236651 

1.38142962 

1.38035113 

1.37913155 

1.37777138 

1.37627114 

1.37463137 

1.37285259 

1.37093537 

1.36888027 

1.36668787 

1.36435877 

1.36189359 

1.35929297 

1.35655756 

1.35368805 

1.35068516 

1.34754960 

1.34428217 

1.34088364 

1.33735487 

1.33369673 

1.32991014 

1.32599604 

1.32195546 

1.31778945 

1.31349913 

1.30908568 

1.30455032 

1.29989437 

1.29511922 

1.29022631 

1.28521720 

1.28009353 

1.27485703 

1.26950954 

1.26405302 

1.25848956 

1.25282137 

1.24705079 

1.24118035 

1.23521271 

1.22915074 


S 

6.00338615 
5.99525141 
5.98589403 
5.97532126 
5.96354074 
5.95056056 
5.93638926 
5.92103581 
5.90450963 
5.88682058 
5.86797902 
5.84799574 
5.82688204 
5.80464968 
5.78131092 
5.75887853 
5.73136578 
5.70478647 
5.67715493 
5.64848603 
5.61879520 
5.58809843 
5.55641230 
5.52375399 
5.49014127 
5.45559255 
5.42012689 
5.38376400 
5.34652429 
5.30842884 
5.26949949 
5.22975880 
5.18923013 
5.14793762 
5.10590624 
5.06316184 
5.01973114 
4.97564182 
4. 930922 50 
4.88560284 
4.83971353 
4.79328640 
4.74635442 


DR/DX 

—.00441553 

-.00520994 

-.00599160 

-.00677546 

-.00755650 

-.00833465 

-.00910987 

-.00988208 

-.01065173 

-. 01141 T 23 

-.01217999 

-.01293943 

-.01369544 

-.01444791 

-.01519671 

-.01594171 

-.01668277 

-.01741974 

-.01815243 

-.01888068 

-.01960427 

-.02032301 

-.02103665 

-.02174496 

-.02244767 

-.02314449 

-.02383511 

-.02451921 

-.02519642 

-.02586637 

-.02652865 

-.02718281 

-.02782838 

-.02846485 

-.02909168 

-.02970826 

-.03031397 

-.03090812 

-.03148997 

-.03205874 

-.03261357 

-.03315354 

-.03367765 


DS/DX 

-.03835181 

-.04517772 

-.05196512 

-.05871169 

-.06541504 

-.07207279 

-.07868250 

-.08524174 

-.09174799 

-.09819875 

-.10459146 

-.11092352 

-.11719229 

-.12339508 

-.12952918 

-.13559180 

-.14156011 

-.14749123 

-.15332222 

-.15907006 

-.16473173 

-.17030404 

-.17578379 

-.18116770 

-.18645239 

-.19163440 

-.19671018 

-.20167607 

-.20652831 

-.21126302 

-.21587621 

-.22036374 

-.22472137 

-.22894467 

-.23302906 

-.23696982 

-.24076200 

-.24440047 

-.24707991 

-.25119473 

-.25433910 

-.25730694 

-.26009184 


CP 

-.01774974 

-.01797955 

-.01820446 

-.01842569 

-.01864130 

-.01885201 

-.01905828 

-.01925893 

-.01945433 

-.01964479 

-.01982960 

-.02000878 

-.02018279 

-.02035090 

-.02051321 

-.02066953 

-.02081977 

-.02096369 

-.02110102 

-.02123204 

-.02135618 

-.02147325 

-.02158351 

-.62168643 

-.02178134 

-.02186890 

-.02194834 

-.02201901 

-.02208175 

-.02213557 

-.02217964 

-.02221494 

-.02223962 

-.02225487 

-.02225953 

-.02225226 

-.02223468 

-.02220515 

-.02216227 

-.02210762 

-.02203949 

-.02195619 

-.02185930 



ORIGINAL PAGE IS 
OF POOR QUALM 


X 

R 

S 

DR/DX 

DS/DX 

CP 

29.70000000 

1.22299747 

4.69895180 

-.03418482 

-.26268709 

-.02174585 

29.88000000 

1.21675617 

4.65111401 

-.03467390 

-.26508562 

-.02161759 

30.06000000 

1.21043031 

4.60287790 

-.03514364 

-.26727996 

-.02147246 

30.24000000 

1.20402363 

4.55420175 

-.03559267 

-.26926223 

-.02130822 

30.42000000 

1.19754012 

4.50536536 

-.03601952 

-.27102408 

-.02112679 

30.60000000 

1.19098405 

4.45617013 

-.03642260 

-.27255663 

-.02092448 

20.78000000 

1.18436002 

4.40673921 

-.03680017 

-.27385044 

-.02070287 

30.96000000 

1.17767296 

4.35711755 

-.03715035 

-.27489541 

-.02045944 

31. 14000000 

1.17092816 

4.30735210 

-.03747109 

-.27568076 

-.02019168 

31.32000000 

1.16413133 

4.25749189 

-.03776015 

-.27619486 

-.01990096 

31. 50000000 

1.15728861 

4.20758825 

-.03601510 

-.27642523 

-.01958318 

31.68000000 

1.15040663 

4.15769496 

-.03823326 

-.27635832 

-.01923943 

31.86000000 

1.14349252 

4. 10786847 

-.03841171 

-.27597947 

-.01886678 

32.04000000 

1.13655402 

4*05816815 

-.03654723 

-.27527266 

-.01846167 

32.22000000 

1.12959949 

4.00865654 

-.03863678 

-.27422036 

-.01802505 

32.40000000 

1.12263800 

3.95939970 

-.03867497 

-. 27260326 

-.01755149 

32.56000000 

1.11567938 

3.91046758 

-.03665895 

-.27100002 

-.01704184 

32.76000000 

1.10873437 

3.86193441 

-.03858342 

-.26878688 

-.01649151 

32.94000000 

1.10181463 

3.81387928 

-.03844299 

-.26613721 

-.01589577 

33*12000000 

1.09493293 

3.76638672 

-.03823165 

-.26302100 

-.01525416 

33.30000000 

1.00810327 

3.71954745 

-.03794256 

-.25940410 

-.01455936 

33.48000000 

1.08134102 

3.67345927 

-.03756805 

-.25524735 

-.01381016 

33.66000000 

1.07466316 

3.62822820 

-.03709922 

-.25050536 

-.01299929 

33.84000000 

1.06808851 

3.58396992 

-.03652566 

-.24512495 

-.01211873 

34.02000000 

1.06163802 

3*54081147 

-.03583600 

-.23904294 

-.01116453 

34.20000000 

1.05533523 

3.49889370 

-.03501551 

-.23216316 

-.01012395 

34.38000000 

1.04920671 

3.45837430 

-.03404732 

-.22445222 

-.00898956 

34.56000000 

1.04328283 

3.41943215 

-.03291049 

-.21573283 

-.00774346 

34.74000000 

1.03759866 

3.38227330 

-.031579 59 

-.20537416 

-.00637143 

34.92000000 

1.03219559 

3.34713997 

-.03001717 

-.19467565 

-.00484739 

35*10000000 

1.02712320 

3.31432446 

-.02817949 

-.16185931 

-.00313552 

35.28000000 

1.02244361 

3.284^9245 

-.02599622 

-.16701783 

-.00118966 

35.46000000 

1 . 31823729 

3.25722582 

-.02336642 

-.14950585 

.00107462 

35.64000000 

1,01461837 

3.23411392 

-.02010512 

-.12817082 

.00379953 

35.82000000 

1.01177372 

3.21600463 

-.01580356 

-.10046594 

,00731804 

36*00000000 

1.01019331 

3.20596552 

-.00878009 

-.05572924 

.01304281 

MACH NO. « 

2.500000 CD 

WAVE - .028562 

SREF * 

6.02628 
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